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Abstract:  This  paper  presents  a  new  approach  for  the  global  elec¬ 
tromagnetic  analysis  of  the  three- Terminal  active  linear  and  nonlinear 
microwave  circuits  using  the  Finite-Difference  Time  Domain  (FDTD) 
Method.  Here,  we  have  updated  the  both  electric  field  components 
on  the  three  -  terminal  active  device  by  correlating  the  voltage  and 
current  with  its  impedance.  This  approach  is  applied  to  the  analysis 
of  a  linear  amplifier  which  includes  a  three-terminal  active  MESFET 
device.  Simulations  results  are  in  good  agreement  with  those  of  the 
commercial  tool. 
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1  Introduction 

The  design  of  today  modern  high  frequency  electronic  systems  requires  the 
inclusion  of  electromagnetic  coupling  effects.  This  can  be  achieved  by  using 
a  full  wave  approach,  which  can  include  the  interaction  between  the  active 
devices  and  their  surrounding  EM  field. 

The  finite-difference  time  domain  (FDTD)  method  [1]  is  one  of  the  best 
candidates  to  be  used  in  high  frequency  system  simulation  model  to  handle 
the  field  effect  in  the  system.  Several  research  works  have  been  reported 
concerning  the  extension  of  the  FDTD  method  to  incorporate  lumped  mi¬ 
crowave  devices  into  3-D  full  wave  analysis.  Among  these,  an  extended  FDTD 
method  [2,  3]  has  been  proposed  to  include  simple  lumped  elements.  How¬ 
ever,  this  approach  may  become  quite  cumbersome  in  the  modeling  of  three- 
ternrinal  active  devices  with  complicated  equivalent  circuit  models  and  seems 
impractical  when  a  multi-port  lumped  circuit  is  involved.  Moreover,  the  ac¬ 
tual  device  size  cannot  be  easily  taken  into  account,  and  tedious  distributing 
procedures  may  be  required  [4].  In  an  effort  to  alleviate  these  problems, 
Thomas  et  al.  [5]  established  the  links  between  the  FDTD  method  and  the 
SPICE  simulator.  Although  this  approach  is  efficient  in  the  sense  that  it 
exploits  the  already  available  codes,  programming  the  link  between  FDTD 
and  SPICE  can  be  time  consuming.  Later  on,  C.  N.  Kuo  et  al.  proposed 
in  [6]  another  approach  where  the  device- wave  interaction  is  modeled  through 
equivalent  voltage-sources  and  state  equations.  In  performing  the  calculation 
of  state  equations,  however,  the  difference  equations  were  not  created  at  the 
same  FDTD  time  step.  Another  approach  has  been  proposed  in  [7,  8]  incor¬ 
porating  a  general  two-terminals  devices  in  a  single  FDTD  cell.  To  the  best 
of  our  knowledge,  this  approach  has  not  been  extended  yet  to  three-terminal 
active  devices.  Thus,  it  is  desirable  that  a  more  straightforward  and  sim¬ 
pler  method  is  developed  to  include  three-terminal  active  devices  into  FDTD 
algorithm. 

In  this  paper,  a  new  approach  is  presented  without  using  state  equations 
which  must  be  solved  at  each  time  step.  Therefore,  it  provides  important 
savings  in  computational  cost  over  [6].  Furthermore,  the  three-terminal  active 
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device  can  be  easily  extended  over  multiple  FDTD  cells  without  a  significant 
modification  to  the  method  and  without  an  increase  in  computational  cost. 

2  Theory 

In  order  to  include  the  equivalent  circuit  model  of  the  MESFET  Fig.  1  into  the 
FDTD  mesh,  the  active  device  can  be  inserted  as  shown  in  Fig.  2  (a),  where 
V\  and  V2  are  aligned  to  the  FDTD  grid  edges  beneath  and  perpendicular 
to  the  microstrip  line  with  each  ends  connected  to  the  microstrip  line  or  a 
grounded  via.  Those  vias  are  modeled  as  perfect  conductors  and  provide  a 
voltage  reference  to  V\  and  V2. 

Assume  that  the  V±  and  V>  are  connected  at  the  nodes  rn  and  vr 

Z1  z2 

respectively  along  the  z-direction;  the  update  equations  for  E™+1(tez1)  and 
E™+1(teZ2)  at  the  two  ports  of  a  MESFET  can  be  found  as: 

Enz+\rEzi)  =  E?(rEzi)  +  ^  [V  x  tf]?+1/2  (rEzi ) 


Fig.  1.  Equivalent  circuit  for  the  JS8851-AS  MESFET 
without  parasitic  elements  (after  [9].  Cgd  = 
0.06  pF,  Cds  =  0.26  pF,  Cgs  =  0.69  pF,  R;  = 
1.42 'D,  gm  =  65  rnS,  gds  =  Rds  =  197 'D. 


microstrip  microstrip 


- ^ 

r“i 

)V2  Z, 

- ►  Y 

Ground  Plane 


(a) 
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(b) 

Fig.  2.  (a)  Implementation  of  the  three-terminal  active 
device.  (b)  structure  of  an  amplifier  without 
matching  network 
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At 

2 7 

nz{rEz 
At  - 
2e  - 


rn+ 1 


m+ 1 


(rEzl  )  +  Jpz(rEzl  ) 

(1) 

—  [V  X  tf]?+1/2 
e 

(r^2) 

(rEZ2  )  +  JFz{rEZ2  ) 

(2) 

The  admittance  matrix  of  a  MESFET  in  frequency  domain,  written  as, 


D(u) 


ju)Ai  —  bj1  A2  juAs  +  lu2A2 

A4  +  jojA5  +  lu2A2  A§  +  jtoA'i  —  u>2Ag 


(3) 


A\  —  CgS  ~t~  Cgd)  ^2  —  R-iC gsCgdj  A3  —  Cgd,  2I4  —  5m 
^5  =  QmRiC gs  —  Cgd,  A§  =  gdsi  Ay  =  Cds  T  Cgd  T  SdsRCgs 
A%  =  RiCgS{Cds  +  Cgd),  D(s )  =  1  +  jUlRiCgs 
The  resulting  admittance  matrix  representation  in  the  this  domain  is  given 

by 

2 

h(u)  =  Y,Yiq(u>)Vq(u)  (4) 

q=l 

2 

h(u;)-Ids  =  J2Y2q(u)vqH  (5) 

q=l 

Where  q  =  1,2  is  the  port  indexe.  The  coefficients  C^q  and  dr  are  obtained 
as  in  [8]. 

By  replacing  jui  with  and  —a ;2  with  J^-,  the  finite  difference  form  of 
Eq.  (4)  is  easily  obtained  as  follow: 

2  2 

d0I?+1  =  J2  C}1VTr+1+  X  °r2v 2n~r+1  ~  diJi 


r= 0 


r=0 


d(iII2l+l  =  Y  C21v1n~r+1+  Y  C22v?~r+1  +  Y  drind7r  -  dill 


(6) 

(7) 


r=0  r=0 

Using  the  following  relationships: 


r=0 


v7  =  -e:{tEzi)Az 
V?  =  -E:(rEz2)Az 

d\  =  -Jiz{rEzl)^xAy 
d2  =  -J^rEjAxAy 

Eq.  (6)  and  Eq.  (7)  become 

dnJ?z+\rEzi)AxAy  =  C^AzE^H^ )  +  C^AzE^HteJ  -  a 
d*J?z+\rEz2)AxAy  =  C21  AzE^1  {r  Ezi )  +  C22AzEnz+\rEJ  -  f?z 
Where  a ^  and  /^represent  the  memory  of  the  circuit: 

a?  =  X  -C^AzE^KteY  -  Y  Cl2AzEn2-r+\rEz2) 


(8) 

(9) 

(10) 


r= 1  r= 1 

-  d\AxAy.Jl,z{rEzi) 

f%  =  X  -cfA^K)  -  x  a22A^rr+i(rEz2)  +  x  4/; 

r= 1  r=l  r=0 

-  diAxAyJ?z(rEz2 ) 


n—r 

ds 


(ii) 
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Substituting  Eq.  (9)  into  Eq.  (1),  Eq.  (10)  into  Eq.  (2)  respectively  and  using 
the  following  definitions: 

_  AtC^Az  _  __AtCpAz_ 

^Z1  2edoAxAy  ’  ^ 21  2edoAxAyqzl 

<*>&,)  =  2-(E?(rE,t)  +  ^  [V  X  ffC+1/2  (raJ  -  ^JK(r&1) 

42:1  \  £  Zt 

A  to" 

H - A  2a 

2ed()AxAy 

_1  AtC^Az  AtC^Az 

^2  2edoAj;Ay’  ^Z2  2edoAxAyqZl 

S^teJ  =  2-fe(rE,J  +  ^|Vx  H]?+1/2  (raj  -  ^J£(  r&J 
2ed()AxAy 

The  discretized  Ampere-Maxwell’s  Eq.  (1)  and  Eq.  (2)  become 

£?+1(r&J  =  7..E?+1(rE,J  +  4”(raJ 
E?+1(r&J  =7«£r‘(rE,)  +  4”(rE.J  (14) 

From  (14),  we  obtain  easily  the  update  equations  for  E™+1(tez  )  and 

£?+1(rx2) 

®?+1(rB,J  =  - 1 - (ij(r&1)  +7xJ?(raJ) 

1  -7^722 

E”+1(raJ  =  1  fe"(raj  +  7a«?(raJ)  (15) 

1  7zi  1Z2 

3  Numerical  results 

We  have  tested  the  validity  of  the  new  approach  on  the  linear  amplifier  con¬ 
figuration  shown  in  Fig.  2  (b),  where  Ids  =  fjmYi  ■  We  selected  the  JS8851-AS 
linear  MESFET  depicted  in  Fig.  1  (a).  To  get  a  50 17  characteristic  impedance 
using  a  substrate  material  of  a  thickness  of  0.795  mm  and  a  relative  dielectric 
permittivity  of  2.2,  the  width  of  the  microstrip  line  becomes  2.4 mm.  The 
cell  sizes  are  Ax  =  Ay  =  0.4  mm,  Az  =  0.265  mm  and  At  =  OAps.  The 
unsplit-field  PML  (crmax  =  60,  m  =  4,  d  =  10)  is  used  to  truncate  the  cir¬ 
cuit.  The  active  region  occupies  four  cells  in  the  y-direction.  Using  this  set 
of  parameters  and  a  Gaussian  pulse  excitation  with  15  GHz  bandwidth,  we 
performed  the  FDTD  simulation. 

Fig.  3(a)  and  3(b)  show  the  magnitudes  of  Sll  and  S21,  respectively, 
simulated  by  our  approach  and  the  Hewlett-Packard’s  Advanced  Microwave 
Design  System  (ADS).  Very  good  agreement  can  be  seen  between  the  two 
sets  of  data  proving  the  validity  of  our  model.  Because  of  the  numerical 
dispersion  of  the  FDTD  simulation  and  also  since  the  vias  were  modeled  as 
perfect  conductors,  there  is  a  small  discrepancy  between  the  two  results. 

4  Conclusion: 

A  novel  approach  is  proposed  to  include  three-terminal  active  devices  into 
FDTD  without  the  need  to  use  the  state  equations  approach.  This  allows 
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(a) 


(b) 

Fig.  3.  The  calculated  results  using  the  proposed  ap¬ 
proach  and  a  commercial  tool  S-Parameters  (a) 
|Sn|.  (b)  \S2i\. 


considerable  savings  on  computational  efforts.  Also,  this  approach  allows  the 
active  device  to  be  easily  extended  across  multiple  cells  without  a  significant 
modification  to  the  algorithm  and  without  an  increase  in  computational  cost. 
Simulation  results  show  very  good  agreement  with  ADS.  Furthermore,  this 
approach  can  be  easily  used  for  the  modeling  of  nonlinear  active  devices. 
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